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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

THEORETICAL INVESTIGATION OF A PROPORTIONAL -PLUS -FLICKER 

AUTOMATIC PILOT 


By Ernest C. Seaberg 


SUMMARY 


The proportional -plus-flicker automatic pilot operates by a non- 
linear principle whereby a fast-acting flicker servomotor response is 
combined with a low-speed proportional servomotor response for the pur- 
pose of obtaining supersonic stability and control. Essentially, the 
autopilot maintains a zero reference about which the output is propor- 
tional to the input. However, a flicker response overrides this propor- 
tional response at a fixed angle of gimbal displacement on either side of 
the zero gyroscope reference. Therefore, in contrast to other high- 
speed control systems, the design requirements are simplified because 
the two components of the proportional-flicker control system are easy 
to build separately and they can be combined in a relatively simple 
manner. 

By application of the proportional-flicker principle, satisfactory 
stability can be obtained by the proper adjustment of the variable fac- 
tors in the autopilot mechanism; namely, the proportional gain, the 
amplitude of flicker control deflection, the autopilot time-lag factor 
(the time lag between flicker and proportional operation), and the point 
in the range that the autopilot switches from a flicker to a proportional 
system. There is a possibility that these factors can be adjusted so 
that a more rapid response time (the time to reach steady state) is 
obtained with the nonlinear proportional -flicker autopilot than with a 
purely linear proportional autopilot. 

For the main part of this analysis, the proportional part of the 
system is approximated by a zero -phase -lag proportional autopilot with 
the assumption that the control surface moves instantaneously at the 
point where the system switches from flicker to proportional. Good cor- 
relation is shown between the results obtained by this method and results 
obtained by using a close approximation of an actual autopilot transfer 
function for proportional autopilot operation. 
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The. proportional -flicker control system, appears to he a practical 
method for obtaining pitch stabilization of a supersonic pilotless air- 
craft. Therefore, trials of this system, particularly in supersonic 
vehicles, appear warranted. 


INTRODUCTION 


As part of the general research program for investigating various 
means of automatic stabilization, the Pilotless Aircraft Research 
Division of the Langley Aeronautical Laboratory has been conducting a 
theoretical analysis to determine the feasibility of using a proportional 
plus-flicker automatic pilot for stabilization and control of a super- 
sonic canard airframe. The principle of operation of this autopilot is 
believed to be unique because it combines a fast-acting flicker servo 
response with a low-speed proportional servo response in a relatively _ 
simple manner, that is, by overriding the proportional part of the 
system at a fixed angle of gimbal displacement from the zero gyroscope 
reference through the use of simple electrical pickoffs attached to the 
displacement gyroscope outer gimbal. 

The reason for attempting to develop an autopilot of this type is 
to overcome the apparent difficulty in building a high-speed proportional 
servomotor by the construction of a proportional -flicker, servo. Because 
the fast-acting flicker servo and the low-gain, slow proportional servo 
by themselves have already been tried and proven, the use of a servo- 
motor combining the two characteristics Is suggested, the fast-acting 
flicker portion to alleviate quickly initial disturbances an'd to secure” 
a large response quickly to an input signal. The amplitude of flicker 
control deflection is larger than available from the proportional part, 
the main function of which Is to secure stabilization around the neutral 
point of the range. 

The analysis contained herein pertains to one specific supersonic 
model configuration for which satisfactory stability was achieved. It 
is very probable, however, that the optimum adjustments of the variable 
factors in the autopilot mechanism have not been realized. It is 
believed that a more complete analysis utilizing an analog computing 
machine would show the optimum autopilot adjustment more closely. The 
results of the analysis contained herein show the effects of the fol- 
lowing conditions on the stability of the autopilot -model combination 
based mainly on the flight condition anticipated as a result of previous 
flight tests of the model: 

1 . Normal acceleration 
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2. Aero dynamic out-of-trim moment 

3* Static margin, altitude, and Mach number variation 

The analysis was continued further to include the response to an 
initial disturbance and the response to a command signal using approxi- 
mate physical autopilot transfer functions in the band of proportional, 
autopilot operation. 


SYMBOLS 


K 


canard control- surface deflection, trailing edge down 
denotes positive deflection* degrees 

error angle between gyroscope frame and instantaneous 
airframe attitude, degrees 


proportional autopilot gain factor 


(*-f) 


e,e 0 


a 

a 

7 

6 


pitch angle measured from horizontal, positive when nose 
is above the horizontal reference, degrees 

first derivative of Q with respect to time, degrees per 
second (d 0 /dt) 

angle of attack, positive when nose is up relative to 
flight path, degrees 

first derivative of , a with respect to time, degrees per 
se cond ( da/dt ) 

flight -path angle, degrees (7 = G - a) 

acceleration due to gravity, 32.2 feet per second per 
second 


n 

M 

Y 


normal acceleration in g units 
Mach number 

stability axis which passes through center of gravity and 
is perpendicular to vertical plane of symmetry 


IY 


moment of inertia about Y-axis, slug feet square 



k 
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L 

m 

4 

5 

c,MAC 

SM 

V 

CL 

Cm 

CL5 

C La 

Cats 

°Ba, 

Cmg 

Cm a 


lift, pounds 

mass, slugsj or pitching moment, foot pounds 

dynamic pressure, pounds per square foot 

wing area, square feet 

mean aerodynamic chord, feet 

static margin, feet 

velocity, feet per second 

lift coefficient (L/qS) 

pitching -moment coefficient (m/qSc) 

SCl/^5 

dCL/dci 

dCm/dS 

dCm/dcc. 

SCm/de 

aCm/Sa 


&t 

m t 

Cmt 

t 


value of control -surf ace deflection which counterbalances 
out-of-trim moment, degrees 

out-of-trim moment caused by model misalinement, foot 
pounds 

out -of -trim-moment coefficient (mt/qSc) 
time, seconds 


t time-lag factor, seconds 

tQ°,t^°,t 1 Q° time that the model attitude crosses the point in range 

denoted by subscript, seconds (for example, to° signi- 
fies the instant of time that the model attitude 9 
is 0°) 
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CD 

c.p. 
c • § • 
D 
s 


L _1 [f(s)] 

Subscripts: 

i 


frequency, radians per second 
model center of pressure 
model center of gravity 
differential operator ( d/dt ) 

Laplace tr ans form variable corresponding to differential 
operator 

signifies that inverse Laplace transform of function f(s) 
is to be taken 


input or forcing function corresponding to a command 
calling for a change in attitude or to a sinusoidal 
input variation 

output or response function, for example, system resporise 
to a command signal or to a sinusoidal input variation 


AIRFRAME AND AUTOPILOT PRINCIPLES AND DESCRIPTION 


A typical qualitative curve of the static servomotor displace- 
ment S plotted against the error signal e for a proportional -flicker 
autopilot system might be as follows: 



<MBB Mil HI I III 
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For the theoretical analysis, the estimated pitch derivatives of 
the canard configuration reported in reference 1 were used. A photo- 
graph and plan-view sketch of the model configuration are shown in 
figure 1. 

The autopilot system, as described in the following paragraph, 
constitutes one method of obtaining a proportional -flicker response, of 
which many variations are possible. It is also feasible that a com- 
pletely electric autopilot system could.be devised which would function 
on the proportional-flicker principle. 

A schematic diagram of the type of automatic pilot investigated in 
this analysis is shown in figure 2. This system consists of a displace- 
ment gyroscope and a rate gyroscope which transmit error signals to a 
diaphragm by means of pneumatic Askania pickoffs. (The signal is trans- 
mitted by an air jet to either of two holes in the pickoff 'block, which 
is connected to the servomotor diaphragm by rubber tubes.) The dia- 
phragm actuates the servomotor slide valve and the servomotor response ‘ 
becomes proportional to the input through the use of a feedback spring. 
The system thus described constitutes the proportional part of the auto- 
pilot, which secures stabilization about the zero gyroscope reference 
point. The flicker portion of the system is obtained by use of elec- 
trical override pickoffs which are mounted on the displacement gyroscope 
frame and which make contact at a preset angle of displacement with a 
pickoff attached to the outer gimbal of the displacement gyroscope. 

(The angle of displacement at which these pickoffs are. set determines 
the switching point arid the width of the band of proportional operation 
about the zero gyroscope reference point.) When either of the flicker 
pickoffs makes contact, one of the override solenoids is energized which, 
in turn, actuates the servo slide valve to cause the flicker action of 
the servomotor piston. The function of the leaf springs in the servo 
feedback linkage is to alleviate the feedback spring and diaphragm 
forces during flicker operation. With this arrangement, more rapid 
flicker action is assured with relatively small solenoids. In operation, 
a time lag occurs at the autopilot switching point. This causes a delay 
between proportional and flicker operation when the airframe attitude 
passes out of the proportional band and causes an overshoot under 
flicker operation when the attitude crosses the band from the other 
direction. It is possible to have delay and overshoot periods of 
different magnitudes. 


METHOD OF ANALYSIS 


The analysis contained herein consists of calculating the tran- 
sient responses of the missile-autopilot combination for various initi al 
conditions, approximating those which might be encountered in flight. 



NACA EM L50I19 


7 


The principle on which the proportional-flicker autopilot operates is 
nonlinearj however, it does lend itself readily to analysis by means of 
linear approximations and a step-by-step solution, the first step of 
which is the response to a constant control deflection corresponding to 
flicker servo action caused by an initial attitude disturbance outside 
of the proportional band. The airframe motion and corresponding auto- 
pilot response preceding the initial disturbance were not considered in 
this analysis. In the initial phase of the calculations, the autopilot 
operation in the proportional band was approximated with a gyro-actuated 

control of the type used for roll stabilization in reference 1. 

Calculations based on this assumption are presented in the results as 
approximations of the actual, autopilot -model pitch transient responses 
to an initial disturbance. Since the actual proportional -plus -flicker 
autopilot cannot incorporate the use of a gyro-actuated control in the 
proportional band because the intricacies of the system necessitate the 
use of a servomotor, further calculations were made using an approximate 

function — - f(D) for the autopilot response in the proportional band. 

Transient responses incorporating this approximation are also presented 
in the results. 


Three methods of analysis were considered in obtaining the pitch 
transient responses of the supersonic model-autopilot combination pre- 
sented in the results. The solution for a constant control- surface 
deflection applies in the first step of each method. The first method 

can be generalized as the response to an initial disturbance (— = K in 

\ 

proportional band), the second as the response to an initial disturbance 
(— = f(D) in proportional band), and the third as the response to com- 
mand signal ( — = f(D) in proportional band). 


Response to an Initial Disturbance 

fb \ 

( — =■ K in Proportional Bandj 


The forms of the equations of motion for constant speed, level 
flight, are: 


5T3i* D - c ^ d 


M 


Cn^ + d) a = CmR 5 


-m5 


mV 

57.3qS 


D 9 - 


mV 

57.3qS 


D + CLqJcc. = Cl S 8 
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The general solution of these equations is broken down Into two steps: 

(l) B = constant for the flicker portion of the solution, and (2) B = K€ 
for the proportional part of the solution. For both steps the methods 
of Laplace (references 2 and 3) were used in obtaining the equation 

of 8 as a function of time [e = f(t)J. 

A qualitative example of a typical solution is" as follows: 



band. Under this condition the control surface B is against its stop 
and the solution for B = constant applies until 9 crosses the pro- 
portional band. At this point a time-lag factor r is introduced to 
allow for the delays in the override solenoid and servomotor. After 
this overshoot period the control surface is assumed to move instantane- 
ously to a value which corresponds to 9 by the relation B = Ke =* -K9. 
The control surface would actually move as shown by the dotted line. 
However, the flicker response is estimated to be at a rate of at least 
700° per secondj therefore, the step approximation is used because the 
solution is simplified, as will be shown, without introducing an appre- 
ciable error. At the end of the delay period t , the proporti onal 
control-surface motion applies and the response is calculated according 
to the relation B = Ke. The initial conditions of this second step are 
obtained from the end conditions of the first step. Since the nature of 
the solution requires that the initial value of a be known, the tran- 
sient ja = f(t£] for a constant 8 was also derived in order to 
determine the value of a at the end. of the first step. 
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Response to an Initial Disturbance 
(jr = f(D) in Proportional Band^ 

This method of solution was derived in order to define more com- 
pletely the motion of the 6 transient at the beginning of the second 
step. The solution for the first step is identical with that of the 
first method of analysis and the derivation for the solution of the 
second step is as follows. 

The relations between the airframe and autopilot parameters gov- 
erning the solution of the second step can be drawn in block diagram 
• form as 



where the transfer functions are functions of the differential operator. 

Using the methods of Laplace (references 2 and 3 ) j the equations of 
motion are transformed to the form 

f ]_( s )0 Q ( s ) + f 2 (s)a(s) = CmgSCs) + f 3 (s) (l) 

flj.(s)e 0 (s) + f^CsWs) = Cl 5 8(s) + fg(s) (2) 

where f 3 (s) and f6(s) contain the initial- condition terms. 

Expanding equations (l) 

0 O ( S ) 

It remains to define &(s). 

6 ± = 0 


and (2) and solving for 9 0 (s) yields 

f7(s)8(s) + f g( s ) 

= £9(3) 

By transforming the relation and using 


5 = f(D)e = -f(D)0 o 
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an equation of the form 

5(s) = f 10 (s)e o (s) + f u (s) (4) 

is obtained where f^(s) contains the initial-condition terms 
involving 5 and ' 0 O . * 

Combining equations (3) and (4). and solving for 0 O ( S ) yields 


0 O ( S ) 


f io ( s ) 
•^*13 ( s ) 


(5) 


where 8 ) contains the required initial values of 9 , a, and 5 

and their derivatives. The inverse transform of equation (5) is the 
transient solution for the second step 


9 ( t) = L -1 



By combining equations (3) and (4) and solving for S(s) an 
equation defining 5 as a function of s can be obtained, the inverse 
transform of which is the second-step solution for the 5 transient 
response. 


Response to a Command Signal 
(j = f(D) in Proportional. Bandj 

In an automatically stabilized missile, a command signal can be 
obtained by changing the gyroscope reference point. This command signal 
may be generated, for example, by the radar unit of a target seeker, "by 
an outside source in a guided missile, or by an altitude control.. When 
the command signal is larger than the proportional band," the flicker 
portion of the autopilot under consideration will function, yielding a 
constant control -surface deflection, and the motion will be initiated 
toward the desired new attitude or new neutral point in .the range. An 
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example of the G 0 variation with time for this third method of 
analysis might be as follows: 



As is shown, the 9q transient response is assumed to start from level 
flight and the new neutral point is outside of the original proportional 
band. Initially* the equations of motion are treated in the same manner 
as in the first step of the first and second method except that the 
initial value of 0 O is zero. The transient response (0 O ) ^ 0 a 

constant 5 applies until 0 O crosses the new proportional band. At 
this point the time-lag factor . ( T ) is again introduced to complete 
the first step. 

The method of the derivation for the solution of the second step 
found herein is carried through in general form in references 2, 4, 
and 5- As mentioned previously, the autopilot frequency response is 
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approximated by a transfer frunction of the form — = f (D) where the 

error angle (e) is the angle between the gyroscope reference and the 
instantaneous airframe attitude. The second step, shown as the dotted 
portion of the 0o transient response, is the response to a step 9 ± 
where 9 ± is the angle between the new neutral point and the value of 8 0 
at the end of the first step, a new zero reference being taken for 0 O 
at this point. The equations of motion are handled in the same manner . 
as for the second step of the second method of analysis except that 
the relation governing the definition of. S(s) is -• 

6 = f (D) e = f(D)(0i - 9 0 ) 

therefore, equation (4) becomes 


5(s) = f 10 

and the solution for 0 O ( S ) 

0 O ( S ) = 



takes the form 


^l 4 (s)e 1 (s) + f 12 (s) 


f ll(s) 


( 6 ) 


The inverse transform of equation (6) gives the solution of the second 
step of the transient response to a command signal. 

The three aforementioned methods of analysis describe the general 
theories which were derived to obtain the solutions for the transient 
responses presented in this analysis. Any deviation from these general 
theories will be described in further detail in the results. 


RESULTS AND DISCUSSION 

Respohse to an Initial Disturbance 
= K for Second Step) 


Determination of the time-lag factor r . - The time-lag factor for 
the autopilot previously described is defined as the time that it takes 
for flicker servo action after the flicker-actuating solenoid has been 
energized or deenergized. In order to obtain an approximate value of 
this factor for use in the calculations, a servomotor was constructed 
which operated on the proportional -flicker principle and which would be 
suitable for actuating the canard control surfaces of a model of the 
type which was flight-tested in reference 1. The results of an 
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experimental test of this servomotor are presented in figure 3* As is 
shown, an average r of 0,03 second was obtained. It is not believed 
that more rapid flicker actuation can be obtained with this type of 
system due to the limitations on the size of solenoids j however, it 
would not be difficult to increase the' time-lag factor if necessary for 
stability. 

Determination of the proportional gain and amplitude of flicker 
control deflection.- For these calculations the coefficients and model 
longitudinal pitch derivatives were estimated on the basis of sea-level 
flight at a Mach number of 1.8 with the static margin equal to 0.86c or 
approximately 14 Inches. The geometric characteristics and estimated 
derivatives used are given in table I. 

A detailed description of the steps required for a solution such 
as shown in figure 4 Is as follows: The first step is the response to 

a constant 6. The second step which starts at some time after the 
model attitude passes the proportional boundary Is governed by the 
approximate relation 5 = -K 9 for the major part of this analysis. 

This simplifies the calculations and still gives a good approximation 
of the actual 9 response. The values of the Inltial^-condltion terms 
necessary for the second-step solution of the equations of motion are 
obtained from their values at the end of the first step. If the 9 
response stays inside the proportional band and the motion dies out, 
the solution is complete in two steps. However, if the model attitude 
again reaches the proportional band, the flicker -actuating solenoid will 
be energized, making more steps necessary. For the third step the 5 
response is assumed to hold a constant value which is determined from 
the relation 6 =--K0 where 9 is the proportional band limit. The 
initial -condition terms are again found from their values at the end of 
the previous step. The time’-lag factor determines the duration of the 
third step. When this period of time has elapsed the fourth step is 
initiated for which the constant flicker 5 applies as in the first 
step. This step-by-step solution is continued until the response of 
the model-autopilot combination is determined. 

The first calculated transient responses, presented in figure 4, 
were based on a flicker 8 of ±10°, a proportional band of ±5° which 
is designated in this and subsequent figures by the long-dashed lines, 
end a time-lag factor of 0.03 second. The initial pitch disturbance for 
the curves presented in this figure was 9° the value of K was -1 
for figures 4(a) and 4(b), while in figures 4(c) and 4(d) K = -0.3 
where K is equal to the control-surface gain ratio 6/e. Figure 4(a) 
shows the response when the first step ends at 0.03 second after the 
model attitude 9 passes the proportional boundary (+5°), designated 
by tej 0 + 0.03 second on the figure. Figure 4(b) shows the response 
when the first step ends at 0.03 second after the model attitude passes 
through the neutral point (0°) or t n ° + ^03 second. This would 
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necessitate a more intricate flicker -actuating pickoff “and relay system 
than that shown in figure 2j however, its development would not propose 
a major problem. As can be seen in both figures 4(a) and 4(b), the 
response is unstable as is indicated by the divergent transient 
oscillations. Curves for the same conditions are shown_in figures 4(c) 
and 4(d) except that K = -0.3 in the proportional band. As can.be 
seen, decreasing the gain in the proportional band does not increase 
the stability but actually has the effect of making the model-autopilot 
combination more unstable. 

Since it is apparent that a flicker control-surface gain of ±10° is 
too high for stability at M = 1.8, calculations for a lower flicker 
gain, namely ±5°, were made. Curves based oh a flicker. S of ±5° and 
on the same derivatives used for figure 4 are. '.presented in figure 5- 
The curves of figure 5(a) show that, for an initial disturbance of 
9 = 9°, proportional band. of ±5°, and for K = -0.3, the” transient 
response is undesirable when the first step ends at t^° + 0.03 second, 
but when the switching point is .extended to to° + 0.03 second the tran- 
sient response is stable. The secondary oscillation wh ic h appears 
during the second step is induced by the instantaneous movement of the 
control surface 8. Figure 5(b) shows that, for a larger K (namely, 
for . K = -1), the proportional band can be smaller (±3°) and the response 
will still be stable and damp out in two steps. For a jailer initial 
disturbance, however, the stability is more critical, as is shown by 
the curve which has an initial disturbance of"±5.5°» For this reason 
the curve shown in figure 5(c) was calculated for an initial disturbance 
which is just outside of the proportional band, namely, an initial dis- 
turbance of ±3.5° for a proportional band of ±3°. As is shown, the cal- 
culated transient response is divergent) however, by increasing the 
time-lag factor t to 0.1 second as in figure 5(^) the transient is 
made to stay inside the proportional band in the second step. Actually, 
the autopilot .in this instance could have two values of r, 0.03 second 
when the attitude passes out of the proportional band, and 0.1 second 
when the attitude crosses 0°. 

Since it appears that a stable transient response can be obtained 
by the proper adjustment of the variable factors in .the autopilot mech- 
anism, it was decided to increase the flicker control -surface gain 
to +7° in order to obtain the advantage of a more rapid response time 
than that obtained with a gain of ±5°. For this reason the analysis of _ 
the calculated transient responses which follows will be based on a 
flicker 8 of ±7° and a proportional band of ±4°. The first responses 
for these conditions, shown in figure 6, were calculated- using K = -1 
in the second step with the first step ending. at some time after the , 
model attitude has passed through the zero gyroscope reference. In this 
instance, the autopilot will function as a proportional system as long 
as the initial disturbance is within, the ±4° boundary. In the event that 
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the model attitude passes either side of these limits, the flicker 
control surface of ±7° will apply, whichever is corrective. The flicker 
actuating solenoid will be deenergized and proportional control will 
again apply some time after the model attitude passes through the zero 
gyroscope reference. Proportional -flicker transient responses to a 9° 
initial 6 disturbance were calculated with the first step ending at 
t q° +0.03 second and at tQ° +0.1 second, as shown in figure 6. In 
this particular instance, the larger time lag yields a transient 
response which is closer to the zero Q reference. The upper curve on 
figure 6 shows the response of a zero phase lag proportional autopilot 

= K for the entire curve^. This curve serves to illustrate the dif- 
ference between the response obtained from a linear autopilot and the 
responses obtained with the proportional -flicker autopilot. For the 
particular value of autopilot gain used in calculating the response 
obtained with the zero phase lag proportional autopilot, the 
proportional -flicker responses, although more oscillatory, are more 
rapid than the pure proportional autopilot response. The value of K 
used for the zero -phase -lag proportional autopilot response was -0.6. 

The choice of this autopilot constant was based on the method suggested 
in reference 5- This method is illustrated in figure 7 where the 
Nyquist diagram with a K of -O .69 is shown to be tangent to, and not 
greater than, the locus of points required to make the peak amplitude 
ratio of the model -autopilot combination 1 . 3 . 

Acceleration effects .- Normal acceleration and 0, a, and 7 
transient responses for the model-proportional -flicker-autopilot com- 
bination are presented in figure 8 . The Q transient is the curve with 
the first step ending at t q° + 0.1 second, as shown previously in fig- 
ure 6 . As is shown, the trim a is about -2.5° and the steady-state 
normal acceleration is approximately -llg for the first step, with a 
peak acceleration overshoot of about -l 8 . 8 g. 

- An indication of the model aerodynamic loads to be expected under 
flicker -autopilot operation can be obtained from the steady-state varia- 
tion of n/6 with' static margin for a constant 6 input, as shown in 
figure 9 for sea-level flight. This curve is based on the variation of 
the model longitudinal derivatives with center-of-gravity location 
at M = 1.8. Figure 9 shows only the steady-state acceleration; he* 'ever 
the peal£ acceleration is the design factor, therefore acceleration 
curves such as the one shown in figure 8 are necessary to estimate the 
amount of acceleration overshoot. The method employed in deriving the 
relation on which the plot of. n/S against static margin is based is 
given in appendix A. 

The effect of out-of-trim moment .- The value of the out-of-trim 
moment is assumed to be that which will give a certain value of 6 
out of trim (5t) or 
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m t = Cm^StqSc 


Is therefore defined as 


This coefficient has the effect of adding another term to the moment 
equation as follows: 

(57.3qSc 1)2 ‘ Cm 0 D ) 9 " ( Cm a + Cm a D ) a = Cm 5 5 + Cm t 

The solution of this equation, combined with the lift' equation, 
for each step of the proportional-flicker transient responses shows the ■ 
effect of introducing an assumed value of an out-of-trim moment. 

Transient responses, including the effect of an out-of-trim moment, 
are shown in figures 10 and 11. These figures are based .on sea-level 
flight at ■ M q = 1.8 and with SM = 0.86c, and the initial 9 disturb- 
ance was k.5 with K equal to -1 for the second step. It is shown in 
figure 10(a) that, without the out-of-trim effect, the model-autopilot 
transient response is made .to damp out in the second step when the time- 
lag factor is increased to 0.1 -second, that is, when the .first step ends 
at 0.1 second after the model attitude has passed through the neutral 
point (0°). Therefore, the remaining .curves shown in figures 10 and 11 
are based on the value tQ° +0.1 second far the length of the first 
step. Figure 10(b) shows the responses when out-of-trim moment coeffi- 
cients of -0.012 and +0.012 are used. These values of 7 Cmt cause the 
6 transients to trim about -0.57° and +0.57°^ respectively. However, 
the stability is not affected to any great extent" since - the response is 
damped to within approximately 1° of the trim value in 0.6 second in 
each case. More severe values of Cmt* namely, +O.O58 and -O.058, were 
used in figure 11. These give trim 6 values of +2.87° and -2.87°, 
respectively, which are only 1.13° from the proportional boundary 
limits (±4 ). However, although the responses take longer to damp to 
the new trim values, they are still stable and stay inside the propor- 
tional band in the second step. For comparison, calculated transient 
responses including an out-of-trim moment and based on a zero -phase -lag 
proportional autopilot with K = -0.6 are also shown in figure 11. For 
these curves a Cmt ±0.035 was used because, for comparison, it is 
desirable that the trim 9 (±2.. 87°) be the same as for the proportional- 

flicker responses. In figure 11(a) it is shown that, since the 
proportional -flicker response must first pass through 0° and then 
oscillate back to the trim value, the response is not much faster than 
the pure proportional autopilot response. However, in figure 11(b) 
where the out-of-trim moment is negative, the proportional -flicker auto- 
pilot has the more rapid response time . 


The out-of-trim moment coefficient ^Cmt) 

Cmt ” Cm§6t 
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The response to an initial di B~burbari.ee including initial values 
of 9 and a, .- An initial pitching velocity of 50 degrees per 

second and an angl e of attack (<x) of 0.5° were used for the results 
presented In figure 12. These are representative of the values which 
might he encountered along the flight range of the type of model being 
used in this analysis. In figure 12(a) the initial $ was 50 degrees 
per second with zero Ini tial angle of attack, while in figure 12(h) 
initial values of both 9 and a were used and the variation of normal 
acceleration and a with time are also shown. The stability of longi- 
tudinal transient responses is not affected greatly by including the 
effect of these initial conditions. 

The effect of Mach number change . - In order to determine the sta- 
bility characteristics of the model-proportional -flicker-autopilot com- 
bination for a different Mach number, the derivatives for M = 1.4 
were estimated as given in table I for the same center -of -gravity loca- 
tion resulting in a static margin of 0.9c. Results based on these 
derivatives are presented in figure 13 as 9 and a transient responses 
to 9° ariri If. 5° initial pitch disturbances. The 9 transients are com- 
parable, except for Mach number, to the curves with the first step 
ending at tQ° + 0.1 second shown previously in figures 6 and 10(a) and, 
although the response time is a little slower, the results of figure 13 
show stability. 

The effect of static-margin reduction . - The constants and aero- 
dynamic derivatives used for the curve presented in figure 14 were based 
on sea -level flight at M = 1.8 and SM = 0.3c and are given in 
table I. This amounts to a decrease in static margin from approximately 
14 inches to 5 inches. As is shown, an undesirable result is obtained 
because the 9 transient response does not die out in the second step. 
Instead, the response diverges until what appears to be a pure flicker 
response is obtained. Therefore, in a system of this type, the problem 
of obtaining a more rapid response time is not solved by simply 
decreasing the static margin. 

The effect of altitude variation .- Except for altitude, the 
proportional -flicker response with the first step ending at to° +0.1 
second, shown in figure 6, is comparable to the curves presented in 
figure 15 where the pitch transient responses to a 9° initial 9 dis- 
turbance are based on flight at 10,000 and 25,000 feet. As is shown, 
flight at altitude produces a slower response time but does not appre- 
ciably affect the stability of the model-autopilot combination because , 
in both cases presented in figure 15, it is indicated that the 
transient response will die out in the second step. 
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Response to an Initial Disturbance 
fj = f(D) in Proportional Bandj 

Determination of the approximate autopilot functions . - The method 
employed to determine mathematical transfer functions which approximate 
the autopilot amplitude and phase response, was as follows: First, the 

experimental autopilot amplitude and phase responses were obtained from 
oscillating-table tests of actual autopilots. Two autopilots were used 
in this analysis. The first consisted of the servomotor used previously 
to obtain the experimental results presented in figure 3 > and the _ 
oscillating-table error angle was measured by a German .V-l displacement- 
and rate gyroscope which generates the input to the servomotor by the 
use of pneumatic Askania pickoffs. The autopilot amplitude and phase 
test points, thus obtained, were then plotted on semilog paper using a 
decibel scale for the amplitude response. Transparent templates based 
on plots of known quadratic and linear functions on this same semilog 
paper were then used until the combination of templates which most 
closely matched the autopilot test points was determined. The method 
used to make the log magnitude and phase templates and explanations of 
their applications are found in chapter 8 of .reference 2. 

The approximate autopilot transfer function obtained by the afore- 
mentioned method is 


5 = 225 (D + 27.2) 

6 (d 2 + 141D + 7TH) 

This function will then govern the proportional pant of the proportional- 
flicker transient responses. A comparison of this approximate function 
with the autopilot test points based on the amplitude and phase response 
of the proportional -flicker servomotor obtained from oscillating-table 
tests of a V-l displacement and rate gyroscope is presented in figure 1 6. 
The agreement between the approximate function and the experimental 
results can be seen in figure 16 (a), where the approximate autopilot 
transfer, function and the autopilot test points, for oscillating-table 
amplitudes of ±1.13°, ±3«03°, and ±4.85° are plotted on„linear coordi- 
nates. The experimental results vary with oscillation amplitude due to 
the nonlinearities of the autopilot mechanism. However, the agreement 
between the mathematical function and the experimental results is sat- 
isfactory except for the low -amplitude oscillations (±1.13°), where the 
amplitude response peaks more sharply. 

A further comparison between the, approximate autopilot transfer 
function and the actual autopilot test points is made in figure 16(b) 
in the form of Nyquist diagrams, where the mathematical function and 
the test points are combined with the model transfer function to obtain 
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the loci of the polar plots. The model transfer function is based on 
sea-level flight at M = 1.8 and with SM = 0.86c. Satisfactory agree- 
ment is again shown between the approximate autopilot transfer function 
and the experimental results except for the Nyquist plot based on the 
±1.13° test points. 

The second approx ima te autopilot mathematical transfer function was 
based on the amplitude and phase response of a servomotor using a dis- 
placement gyroscope only to generate the input si gn a l . Vlhis approximate 
function is 


£ = 2880 
6 D 2 + I56D + 3600 

and it is compared with the autopilot consisting of a displacement 
gyroscope plus servomotor in figure 17. Autopilot test points are 
plotted for oscillating-table amplitudes of ±it-,2 0 and ±2.15 > and satis- 
factory agreement with the plot of the approximate mathematical function 
is shown. 

Autopilot containing displacement -plus -rate gyroscopes . - As shown 
previously in the method of analysis section, the function 8(s) and 
the initial values of the 5 derivatives are included in the derivation 
of the 0 o (s) function in the second method of analysis. Although the 
manual solution for the 9 and S transient responses are much more 
involved when using this type of derivation, it is of value because it 
yields a more complete definition of the control -surface motion. 

Figure 18 shows the pitch transient responses to different initial 
9 disturbances based on the longitudinal derivatives given in table I 
for sea-level flight at M = 1.8 and with SM = 0.86c. The flicker 
gain is ±7° and the proportional band is ±^+° . In figure i8(a), the 
initial 9 disturbance is 10° and the first step ends at 0.03 second 
after the model attitude has passed through the neutral point (0°). In 
this figure the calculated 5 response is also shown. In figure 18(b) 
the initial 9 disturbance is b.^° with curves shown for the first step 
ending at tQ° +0.03 second and to 0 + 0.1 second. The 6 response is 
plotted for the case where the first step ends at to° + 0.03 second. 

As is shown, the s mal ler time-lag factor (0.03 second) is more critical, 
with the second oscillation of the 9 transient actually crossing the 
proportional band. This would ordinarily necessitate handling the 
solution with more than two steps. However, it would be necessary to 
know the initial value of a for the third step and, since the charac- 
teristic equation for the a transient when j = f(D) is tenth order, 

a manu al solution was not attempted. It is more desirable to obtain a 
less os cilla tory transient response which damps out in two steps by 
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increasing the time-lag factor, as is shown for the curve where 
r = 0.1 second. Figure 18(c) shows the 9 response to a 15° initial 
disturbance, where for a -T of 0.03 second the motion is damped in two 
steps . 

The agreement between the results shown for the second method of 
analysis using the autopilot containing displacement-plus-rate gyro- 
scopes and the results shown for the first method of analysis is good. 
The secondary oscillation which appeared in the second step of the first 
method of analysis, for example, is also present when using this second 
method of analysis, resulting in transient responses, the general 
shapes of which are comparable. On this basis the validity of using 

the much simpler relation — = K for the major part of this analysis 
seems to be justified. 

Autopilot containing displacement gyroscope only .- The second step 
of the pitch transient responses to 10^ and 4.5 U initial 9 disturb- 
ances presented in figure 19 is based on the transfer function 

5 = g880 

6 D 2 + 156D + 3600 

These results show that stability can be obtained with a proportional- 
flicker autopilot made up of a servomotor and^di splacement gyroscope 
only. However, a comparison of figures 18 and 19 shows that the auto- 
pilot without rate yields a more oscillatory response in the 
proportional band. 


Response to a Command Signal 
^ ■= f(D) in Proportional Band^ 

The response to a command signal . - Pitch transient responses to 
1.5° and 10° command signals are shown in figure 20. These responses 
are based on sea-level flight at M = 1.8 and with a static margin 
of 0.86c. The initial disturbance in each case is larger than the 
proportional band* therefore, the flicker control setting of +7° applies 
for the first step of the solution, which ends at 0.1 second after the 
transient crosses the desired new attitude. In the second step the 
approximate transfer function 

8 = 225(D + 27.2) 

e (d 2 + lllD +7744) 
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applies. In figure 20(a) the desired change in attitude is 4.5 or, as 
explained in the method of analysis section, the transient response 
shown is the response to a 4.5° command signal. In figure 20(b) , the 
command signal calls for a 10° change in attitude. These results show 
that stable responses to command signals can be obtained with the 
proportional-flicker autopilot and the responses appear to be similar, 
except for direction, to the results for 4.5° and 10° initial disturb- 
ances shown in figure 18. A comparison between the proportional-flicker 
and the zero phase lag proportional response is made for both curves 
shown in figure 20 and in each case the proportional-flicker autopilot 
has the more rapid response time. 


CONCLUSIONS 


The proportional-flicker automatic pilot operates on a nonlinear 
principle, whereby a high-speed flicker servomotor response is combined 
with a low-speed proportional servomotor response for the purpose of 
obtaining stability and control in supersonic flight. Physcially, the 
autopilot motion operates about a zero reference within two bands. In 
the inner band, the autopilot output is proportional to the input and a 
flicker response overrides the proportional response at a fixed angle of 
gimbal displacement on either side of the zero gyroscope reference. The 
conclusions arrived at as a result of the analysis conducted herein, 
based on a specific supersonic model configuration, are as follows: 

1. Satisfactory stability can be obtained by the proper adjustment 
of the variable factors in the autopilot mechanism, namely, the propor- 
tional servo gain, the amplitude of flicker control deflection, the 
autopilot time-lag factor, and the point in the range that the autopilot 
switches from a flicker to a proportional system. 

2. A reasonable aerodynamic out-of-trim moment of the model will not 
affect the stabilization qualities of the proportional-flicker autopilot 
to any great extent. Decreasing the static margin appears to have more 
of an effect on the stability of the proportional-flicker autopilot -model 
combination. For the same autopilot characteristics, a decrease in Btatic 
margin from 0.86 to 0.3 mean aerodynamic chord yields a transient 
response which does not die out but which diverges until what appears 

to be a pure flicker response is obtained. 

3. Good agreement is shown between the response to an initial dis- 
turbance using a mathematical transfer function to approximate the auto- 
pilot and the method of assuming that the autopilot has a pure propor- 
tional response in the proportional band with instantaneous movement of 
the control surface at the switching point. 
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4. The proportional -flicker control system can he fabricated and 
appears to be a practical method for obtaining pitch stabilization of 
a supersonic pilotless aircraft. Therefore, trials of this system, 
particularly in supersonic vehicles, appear warranted. 


Langley Aeronautical Laboratory 

Rational Advisory Committee for Aeronautics 
Langley Air Force Base, Va. 
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APPENDIX A 


DETERMINATION OF THE RELATION WHICH EXPRESSES n/8 
AS A FUNCTION OF STATIC MARGIN 


The relation 


Cm$ = c Ls c 


(x - SM) d£ (y + SM) 

c d 6 c 


(1) 


expresses Cmg as a function of static margin. 

From the lift equation, the steady-state value of the normal 
acceleration in g can he derived, for a constant control-surface 
deflection as 


n 


1 

32.2m. 


-Ci^qS 


F— + c L 5 <lS 
Cm a 


6 


where the steady state a/6 has been assumed to be equal to 
incorporating the relation 



Cm q = SM 
CL a" c 

the following relation is. obtained, where the variation of the steady- 
state n/6 is expressed as a function of SM and Cmg 


n _ qS 
6 32.2m 

For a particular value of SM, Cmg can be evaluated from equa- 
tion (l). Then, by using these two values in equation (2), n/6 is 

determined . * 


Cmg 


SM 


c Ls 


( 2 ) 




!^»g»irxoyi 



2k 
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cl 5c 


CLq W 

x 


7 

de/d& 


Symbols not Previously Defined 


rate of change of canard lift coefficient with canard 
deflection SCl/c)6 c 


lift -curve slope of main wing 




distance from center of pressure of canard control surfaces 
to center of pressure of model, inches 


distance from the center of pressure of main wing to center 
of pressure of model, inches. 

rate of change of downwash angle at wing due to deflection of 
canard control surfaces 


SM static margin, negative when c.p. is behind the c.g.,' inches 
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TABLE I 

ESTIMATED LONGITUDINAL DERIVATIVES 


[All derivatives in degree measure j 
m = 4.922 slugsj c = 1.395 £tj 


Iy =• 37.66 slug -ft 2 > 
s = 2.52 ft 2 [ 


Static 

margin 

Mach 

number 

Cma 


Cmg 

C ma 

CLe 

°La 

0,86c 

1.8 

-0.000138 

-0.052 

0.204 

-0.0000116 

-0.000051 

0.0607 

0.9c 

1.4 

-.000187 

-.0621 

.023 

-.0000205 

-.00105 

.0692 

0.3c 

1.8 

-.0000784 

-.0181 

.0182 

-.0000051 

-.000051 

.0607 


VARIATION OF FLIGHT .CONDITIONS 


Altitude 

(ft) 

Mach 

number 

q 

(lb/ft2) 

V 

(ft/sec) 

sea level 

1.8 

4270 

1963 

sea level 

1.4 

2902 

1562 

10,000 

1.8 

3294 

1937 

25,000 

1.8 

1790 

1835 
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c.p. location 
('station) 

80 42 at M-/.8 
81.03 atM=!4 

(h) Plan -view sketch of model conf iguration. . 

£ 


Figure 1.- Concluded. 


NACA EM L50I19 



30 


NACA RM L50I19 



Figure 2.- Schematic diagram of a proportional- flicker autopilot system. 



Flicker actuating s/gnot Servo displacement , m. 
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Proportional- flicker servo response to a sinusoidal input signal 
obtained from experimental tests. 


Figure 3* - 





j 3 .4 J .6 0 ./ 1 .3 A .S .6 .7 .8 

t , sec fjSec 


(b) of first step at /b) Hud. of first step at 

t^o + 0.03 second, t Q0 + 0.03 second, 

K = -1 for second step. K ~ -1 for second step. 


Figure 4.- Pitch, transient responses to a 9° 

i xr\ x ^4. ~ 4 . \m ^ n Q cm . 

wu aca-j.cvci ixx^uu cvu j*l = x.vj ojjlul uk»i * 

T =s 0.03 second, flicker 5 = ±10°. 
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(d) End of first step at 
tQ° + 0.03 second, 


K = -0.3 for second step 
ciuded. 




- End sf first step at +.03 sec 

- End of first step at t 0 o + .03 sec 


oo 


B\ 












o •/ .2 3 4 S .6 O ./ j? .3 .4 .£ .6 .7 

T , sec t , sec ' 


\.0j .initial 6 disturbance - 3*5''^ t = 0.03 second., ^aj Initial 6 disturbance = 3.5', r = 0.1 second, 
proportional band = ±3°, end of first step at proportional band = ±3°> end of first step at 

t n o + 0.03 second, K = -1 for second step. t rt0 + 0.1 second, K - -1 for second step. 

u 


Figure 5* - Concluded. 
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J 4 .5 .6 .7 .8 .7 /jO /./ A2 /.S 1.4 

T , see 

Figure 6.- Pitch tr ans ient responses to a 9° initial 6 disturbance based on 
sea-level flight at M = 1.8 and SM = 0.86c. Proportional band = ±4°, 
flicker 5 = ±7°, K = -1 for second step. The calculated response of 
'a zero-phase- lag proportional autopilot is show for comparison. 
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Figure 7.- Nyquist diagram for a zero-phase-lag proportional-autopilot model 
combination with K = -O .69 based on sea-level flight of the model and 
the estimated derivatives for M =s 1.8 with SM = 0.86c. 



* 



1 



Figure 8.- Normal acceleration n and longitudinal 6 , a, and 7 transient 
responses to a 9° initial 9 disturbance "based on sea-level flight at 
M = 1.8 and SM = 0.86c. Proportional "band = ±4°, flicker & = ±7°, 
end of first step at t QO + 0.1 second, K = -1 for second step. 
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(b) Comparison of proportional -flicker autopilot with a zero -phase -lag 
proportional autopilot ^ = K j for a negative out -of -trim mcanent. 

Figure 11,- The effect of an out-of-trim moment on the pitch transient 
responses to a 4.5° initial 0 disturbance based on sea-level flight ■ 
at M = 1.8 and SM = 0.86c. Proportional band = ill- 0 , flicker 
B « ±7°j first step ends at tQO + 0.1 second, K = -1 for second 
step . 
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(a) 0 ( 0 ^ = 50 deg/sec, cl( q \ = 0. (b) 0^ = 50 deg/sec, = 0.5°. 


Figure 12.- Longitudinal transient responses to a 4.5° initial 6 disturbance 
including initial values of B and a based on sea- level flight at 
M = 1.8 and SM = 0.86c. Proportional band = flicker 8 = ±7°> 

-P4-«_4- ^4-^, «+ +«r\ 1 ri 1 csonnr^ TT = -1 ppnnnrl 
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0 J 4 3 4 S .6 .7 A .9 /a U U i.3 f4 

T , sec 

(b) Initial 0 disturbance = 4.5°. 

Figure 13.- Longitudinal transient responses to initial 0 disturbances 
of 9° tod 4.5° based on sea-level flight at M = 1.4 and SM = 0.9c. 

1 : Proportional band = ±4°, flicker 8 = ±7°, end of first step at 

tQ 0 +0.1 second, K = -1 for second step. 
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Figure l4. - Pitch, transient response to a 10° initial 9 disturb' 
ance based on sea- level flight at M = 1.8 and SM = 0.3c. 
Proportional band = ±1°, flicker 5 = ± 7 °, end of first step 
at t, o + 0.05 second, K = -1 for second step. 








(a) Amplitude and phase response of the approximate autopilot function 
and of the actual autopilot test points. 


Figure 1 6 .- Comparison of the function -which appro x imates the proportional 
part of the proportional- flicker autopilot with the autopilot test, 
points based on the amplitude and phase response of a proportional- 
flicker servo obtained from oscillating-table tests of a V-l displace- 
ment and rate gyroscope. 
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-270 


—Approximate autopilot 
Autopilot test points 

+ 1 . 13 ° 

— 13.03° 

± 14 . 85 ° 



(b) Nyquist diagrams for approximate autopilot function-model combination 
and for longitudinal pitch oscillation amplitudes of-il.l3°, - 3 * 03 °, 
and ±^.85° 'based on sea-level flight at M = 1.8 and SM = 0.86c. 


Figure 16.- Concluded 
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Figure 1?.- Comparison of the function -which approximates the proportional 
part of the proportional- flicker autopilot with the autopilot test 
points based on the amplitude and phase response of the servo obtained 
from oscillating-table tests of a displacement gyroscope only. 






0 


■ J ■ l - J I 1 — L— L 

J 2 J 


.4 


f- .-I i — . — J, 


-1 1 -J — i-l L I I 1 I I I 1 _L 


.7 £ 

T ,Scc 


.? 


/J 


// /J 


/J 14 


/ V 


(a; initial & disturbance = !U U , ena or nrst step at tgo + U.U3 second. 


toi i O nj j i* j — j l. j — — - J (14 ^ n n jj i i *» j 

J?J.gUX« 0-U* - JTJL UUli OrtUltJXeJU. U X-CHjyUIlHefcl UU tui XJLU- UXOuL U U_L0 XiLUTUtlLLCe uaeeu 

on sea-level flight at M = 1.8 and SM s= 0.86c. Proportional hand = ±4°, 


_ j-v° 8. _ Wn\ 


rate signal. 


• 7 


-P/->iT* OQr>rrt1rl d+nm 4 TirO mH n rt J-l «r\1 n/»iMnAri+_-rtl iia 


. \-** / J.V* UVWW**«. u wj^ 


ku,vXu5 uiDJ^XdiU 




1 is 


NACA RM L50I19 ■ 







NACA RM 




nn.’KiMiCE.v 


«T, 




Figure 19»- Pitch transient responses to 10° -and 4.5° initial 
0 disturbances based an sea-level flight at M = 1.8 with 
SM = 0.86c. Proportional hand = ±4°, flicker B = ±T°, end 

of first step at t^o t Oil second, — = for second 

step using displacement signal only. 
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(t) New neutral point = 10°, end of first step at t^o + 0.1 second. 

Figure 20.- Pitch transient responses to 4.5° and 10° command signals 
based on sea -level flight at M = 1.8 with SM = 0.86c. 

Proportional band = ±4°, flicker 8 « ±7°, ~ = f(D) for second 

step including dlsplacenent-plus-rate signal. The responses of a 
zero-phase-lag proportional autopilot are shown for comparison. 
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